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The investigation focuses on development of a rapid multidisciplinary analysis and optimization capability 
for launch-vehicle design. Two multidisciplinary optimization strategies in which the analyses are integrated in 
different manners are implemented and evaluated for solution of a single-stage-to-orbit launch-vehicle design 
problem. Weights and sizing, propulsion, and trajectory issues are directly addressed in each optimization process. 
Additionally, the need to maintain a consistent vehicle model across the disciplines is discussed. Both solution 
strategies were shown to obtain similar solutions from two different starting points. These solutions suggests that 
a dual-fuel, single-stage-to-orbit vehicle with a dry weight of approximately 1.927 x 10 5 lb, gross liftoff weight of 
2.165 X 10 fi lb, and length of 181 ft is attainable. A comparison of the two approaches demonstrates that treatment 
of disciplinary coupling has a direct effect on optimization convergence and the required computational effort. In 
comparison with the first solution strategy, which Is of the general form typically used within the launch vehicle 
design community at present, the second optimization approach Is shown to be 3-4 times more computationally 
efficient. 


Nomenclature 

A , ~ nozzle exit area, ft 2 

c - nonlinear constraint vector 

F z = wing normal force, lb 

GLOW = vehicle gross liftoff weight, lb 

/ sp = specific impulse, s 

J = objective function 

MR = mass ratio 

= reference aerodynamic surface area, ft 2 
T = thrust, lb 

T/W = thrust-to-weight ratio 

x - design variable vector 

a = angle of attack, deg 

£ = discipline compatibility tolerance 

Subscripts 

c = computed 

vac = vacuum 

si = sea level 

Introduction 

F OR several years, various Earth-to-orbit transportation options 
have been examined with the goal of reducing operating costs 
relative to the current U.S. launch fleet. 1,2 Many of these solutions 
have focused on fully reusable systems employing various levels of 
advanced technology. 3 Although a wide range of options have been 
examined, including single- and two-stage systems using rocket 
and/or air-breathing propulsion, current emphasis has been placed 
on single-stage-to-orbit, rocket-powered vehicles. 3-5 The design of 
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such a vehicle is a multidisciplinary process in which aerodynamics, 
propulsion, weights and sizing, structures, performance, heating, op- 
erations, and cost must be addressed. 5 Although it is imperative that 
each of these disciplines be addressed at the conceptual design level, 
it is equally vital to be able to perform this multidisciplinary analysis 
and optimization rapidly so that the numerous design options may 
be evaluated and understood. 

The present investigation focuses on development of a rapid mul- 
tidisciplinary analysis capability for launch-vehicle design. The spe- 
cific application chosen is that of a dual-fuel, single-stage-to-orbit 
launch vehicle. Two multidisciplinary optimization strategies are 
implemented and evaluated for the solution of this problem, and dif- 
ferences among the approaches are highlighted. Weights and sizing, 
propulsion, and trajectory issues arc directly addressed in the opti- 
mization processes. Additionally, the need to maintain a consistent 
vehicle model across the disciplines is discussed. 

Problem Definition and Disciplinary Analyses 

In this analysis, design of a single-stage-to-orbit launch-vehicle 
includes specification of the ascent trajectory through the initial 
launch azimuth, time of flight, and pitch-angle history. Determina- 
tion of the appropriate component weights and sizes is performed, 
and the vehicle dry weight (defined as the vehicle weight with- 
out payload, propellant, fluids, or crew) is selected as the min- 
imization variable. Table 1 summarizes the design variables and 
constraints that characterize this optimization problem. Propulsion- 
system characteristics to be optimized include the liftoff thrust-to- 
weight ratio, two nozzle area ratios, and two fuel-to-oxidizer mixture 
ratios. In this analysis, a dual-position nozzle is used to provide the 
performance benefit of a smaller nozzle exit area at liftoff (to maxi- 
mize sea-level thrust) while allowing for a larger expansion at high 
altitudes (to maximize vacuum thrust). Two mixture ratios require 
specification because the vehicle is operated in different propul- 
sive modes. Beginning with liftoff, hydrogen and kerosene are both 
burned as fuel, but during a later portion of the ascent, hydrogen 
becomes the sole fuel. Such a dual-fuel strategy has been shown to 
provide significant dry-weight reductions. 6 As listed in Table 1 , the 
transition time from mode 1 to mode 2 propulsion and the extension 
of the dual-position nozz re also optimally determined. 

The vehicle is sized l iver and return a 25,000-lb payload to 
the Space Station folio' aunch from the Eastern Test Range at 
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Table 1 Single-stage-to-orbit launch-vehicle design characteristics for iterative-loop and sequential 
compatibility-constraint methods (objective function = dry weight) 


Iteration 

Design variables 

Min. 

Max. 

Initial value 
1 2 

l 

GLOW," lb 

1 5e6 

3.0e6 

2.0e6 

2.5e6 

2 

f ft2 

2.5e3 

5.0e3 

4.0e3 

3.0e3 

3 

Landed weight,* lb 

1.0e5 

3.0e5 

2.2e5 

2.0e5 

4 

Base diameter,* ft 

20.0 

40.0 

30.0 

25.0 

5 

Launch azimuth, deg 

20.0 

60.0 

45.0 

40.0 

6 

Vehicle liftoff T/W 

1.2 

1.5 

1.3 

1.4 

7 

% LH 2 , mode 1 

3.0 

15.0 

6.0 

4.0 

8 

Mode 2 mixture ratio 

3.0 

15.0 

7.5 

9.0 

9 

Nozzle area ratio, retracted 

40.0 

100.0 

60.0 

50.0 

10 

Nozzle area ratio, extended 

60.0 

200.0 

120.0 

140.0 

11 

Propulsion mode transition Mach number 

2.0 

20.0 

6.0 

5.0 

12 

Nozzle transition Mach number 

0.0 

8.0 

1.0 

2.0 

13 

Fz boundary duration, s 

20.0 

150.0 

80.0 

100.0 

14 

Total trajectory time, s 

300.0 

500.0 

380.0 

400.0 

15-40 

Set of pitch angle, b deg 

-1120.0 

0.0 

Quadratic 

Linear 

Iteration 

Constraints 

Min. 


Max. 

i 

Terminal altitude, ft 

3.03505e5 


3.04105e5 

2 

Terminal velocity, ft/s 

2.5843e4 


2.5845e4 

3 

Terminal flight-path angle, deg 

— l.Oe — 03 


l.Oe - 03 

4 

Terminal inclination, deg 


51.5 


51.7 

5 

Dynamic pressure, lb/ft 2 


0.0 


1000.0 

6 

Fz magnitude, lb 


0.0 

2.5 x 

landed weight 

7 

Nozzle extension ratio 


0.0 


2.0 

8 

a , deg 

- 

-15.0 


15.0 

9 

Pitch rate, deg/s 


-5.0 


-5.0 

10 

GLOW compatibility, 0 lb 

- 

L0e4 


1.0e4 

11 

S rc f compatibility, 0 ft 2 


-5.0 


5.0 

12 

Landed weight compatibility, 0 lb 

- 

1.0e2 


1.0e2 

13 

Base diameter compatibility, 0 ft 


- 0.1 


0.1 


*Not treated as a design variable in the iterative-loop approach (Fig. 4). 
h \n this system, 0 deg is vertical; ±90 deg is horizontal. 

L Not treated as a constraint in the iterative-loop approach (Fig. 4). 


the Kennedy Space Center. For this analysis, the Space Station is 
assumed to be in a 220-n-mile-altitude, circular orbit with a 51.6- 
deg inclination. The single-stage-to-orbit vehicle is initially inserted 
into an orbit with an altitude of 50 x 100 n mile with the correct 
inclination. Onboard propellant is then used to transfer to and circu- 
larize at 220 n mile. This maneuver is accomplished by two engines, 
each with a vacuum thrust of 6000 lb and vacuum 7 sp of 462.2 s. 

Terminal constraints on altitude, velocity, flight-path angle, and 
inclination are enforced, as are maximum in-flight normal-force, 
angle-of-attack, and dynamic-pressure limits (see Table 1). Each of 
these in-flight constraints is included to reflect an analysis restriction 
derived outside the scope of the present investigation. For example, 
a limit on the allowable normal force is included such that the as- 
cent wing loading does not exceed structural limits (2.5 times the 
landed weight) derived in an entry sizing analysis. Furthermore, an 
angle-of-attack constraint is included such that the validity range of 
the aerodynamic analysis is not exceeded. Sizing restrictions require 
that a limit be placed on the allowable extension of the dual-position 
nozzle such that a maximum 2-to-l increase in exit area results. A 
pitch-rate limit of 5 deg/s (to reflect control issues that are not mod- 
eled) and a maximum 3 g acceleration limit are also imposed. Note 
that since the maximum acceleration is enforced within the trajec- 
tory analysis by engine throttling, not by the optimizer, a maximum 
acceleration limit is not listed in Table 1. 

To simplify the analysis so that the problem is tractable, sev- 
eral design disciplines were decoupled from the present analysis. 
An existing vehicle geometry, aerodynamics database, and internal 
packaging analysis were used. 7 Data from aerodynamics, structures, 
heating, and other subsystems were fixed or scaled appropriately. 
Of major significance for the current optimization study are the dis- 
ciplines of propulsion, performance, and weights and sizing, each 
of which is modeled here. 

For this analysis, propulsion-system parametrics supplied by 
Pratt and Whitney, based on modification of the proposed Russian 
RD-701 dual-fuel engine, were used. 8 After a regression analysis, 


‘ Nozzle area ratios 


- Fuel/oxidizer mixture 
ratio, mode 2 


% LH2, mode 1 


T 

For each f T vac 
engine/propellant < Ae 
setti "9 [ lsp va c 



% RP, mode 1 
% LOX, mode 1 
T/W | engine 
Ispsi 

Tsl 

Propellant bulk density 


Fig. 1 Propulsion analysis inputs and outputs. 


this parametric data set is used as shown in Fig. 1. Given the noz- 
zle area ratios and fuel-to-oxidizer mixture ratios, numerous engine 
parameters are computed. These parameters include the sea-level 
engine thrust-to-weight ratio, sea-level specific impulse, and pro- 
pellant bulk density (required inputs to the weights and sizing anal- 
ysis), as well as the vacuum thrust, vacuum specific impulse, and 
nozzle exit area (required trajectory inputs). 

In this investigation, the three-degree-of-freedom equations of 
motion were numerically integrated with use of the Program to Op- 
timize Simulated Trajectories (POST). 9 The vehicle is treated as a 
point mass; Earth rotation and oblateness are modeled, and the 1976 
standard atmosphere is used. As shown in Fig. 2, the required set of 
POST inputs includes vehicle parameters (e g., gross liftoff weight, 
vehicle thrust-to-weight ratio at liftoff, aerodynamic coefficients, 
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r~ GLOW, vehicle T/W | |j ft0 ff, 
landed weight, base diameter 
Initial launch direction, Sref 
Pitch angles 

T vac’ A e- ls Pvac settings 
Propulsion/engine transition Mach numbers 
Fuel/oxidizer mixture ratio, mode 2 
% LH2, % RP, mode 1 

— Other trajectory inputs 

• Aero coefficients 

• Atmospheric model 

• Mission information (orbit, launch site, 
flight duration, normal-force boundary 
duration) 



Fuel/oxidizer 
mass fractions 


Evaluation of _ GLOW 

inflight and Burnout weight 

terminal conditions 

Fig. 2 Trajectory analysis inputs and outputs. 

— 77W | engine^ * s PsJ> ls Pvac 
Tsb T vac 

Vehicle T/W | liftoff 
Fuel/oxidizer mass fractions 
% LH2, % RP, mode 1 
MR 

— Other weights & sizing inputs 

• Vehicle geometry 

• Propellant densities 

• Mission information 
(crew size, duration, orbit) 



Dry weight GLOW c 

Sref c 


Landed weight c 
Base diameter c 

Fig. 3 Weight and sizing analysis inputs and outputs. 

and reference area) as well as trajectory parameters (pitch-angle 
history, launch azimuth, and propulsion-system transition Mach 
numbers) and propulsion parameters (propellant mixture ratios and 
engine performance characteristics). POST is used to evaluate the 
in-flight and terminal constraints and to compute the vehicle mass 
ratio and required propellant fractions (weights and sizing inputs). 

For a given mass ratio, the configuration sizing program (CON- 
SIZ), developed at NASA Langley Research Center, is used to size 
the vehicle and determine the dry weight. Within CONSIZ, the vehi- 
cle is modeled as a collection of components representing structure, 
subsystem, and propulsion elements. Specific weight-estimating re- 
lationships for each component are considered part of the program 
input. This provides user flexibility in both vehicle modeling and 
analysis depth. As shown in Fig. 3, the liftoff thrust-to- weight ra- 
tio, mass ratio, and propellant fractions, as well as several other 
propulsion system parameters, are required inputs to CONSIZ. Note 
that, with the exception of the liftoff thrust-to-weight ratio, all of 
these inputs are computed by one of the other two disciplinary 
analyses. In addition to dry weight, CONSIZ computes the gross 
liftoff weight, reference aerodynamic surface area, base diameter, 
and landed weight (each of which is a required trajectory input). 


From Figs. 1-3, it is clear that solution of this problem requires 
an iterative approach because the trajectory and weights and sizing 
analyses each require inputs that are computed by another disci- 
pline. For example, one must ensure that the reference aerodynamic 
surface area resulting from the vehicle sizing process (5 reff . of Fig. 3) 
is the same as the reference aerodynamic surface area used to com- 
pute the aerodynamic forces and evaluate the flight path (S^, of 
Fig. 2). Consistency must also be maintained in regard to the gross 
liftoff weight, base diameter, landed weight (which effects the as- 
cent normal-force limit), mass ratio, propellant fractions, and each 
of the propulsion discipline outputs. Furthermore, the liftoff thrust- 
to-weight and fuel-to-oxidizer mixture ratios must be treated in a 
similar fashion, since these parameters are input to more than one 
discipline. 

Optimization Approaches 

The disciplinary analysis tools used in this study were originally 
created as independent programs, each operated by a disciplinary 
expert. To obtain a feasible point design vehicle, a design team 
was required to iterate manually among these disciplinary analyses. 
In many cases, “optimization” was performed through trade stud- 
ies in which the parameters were varied one at a time. Significant 
improvement over this one-variable-at-a-time, parametric approach 
has been achieved using response-surface methods (RSMs). 5-7 In 
this multidisciplinary strategy, feasible designs are computed at nu- 
merous statistically selected points in the design space, and a sur- 
face is fitted to these points. Optimization is then performed on 
this approximate representation of the design space. As another 
means of comparison, the launch-vehicle design problem discussed 
in the present investigation is solved using system sensitivity anal- 
ysis (SSA) in Ref. 10. The present results are compared with those 
obtained with SSA in that reference. 

The present investigation considers two solution approaches 
in which the previously independent disciplinary tools are inte- 
grated. These two optimization strategies will be referred to as 
1) the iterative-loop approach and 2) the sequential compatibility- 
constraint method. For this particular problem and set of analysis 
tools, these integrated approaches provide several advantages over 
the use of RSMs. These advantages include the following: 1) op- 
timization over the actual design space as opposed to the use of 
an approximate surface, 2) no human interface time requirement to 
negotiate between disciplinary programs, and 3) the guarantee of a 
consistent multidisciplinary vehicle model at the solution. The pri- 
mary disadvantage to the use of either of these integrated strategics 
is the up-front time required to integrate the disciplinary tools in a 
general manner. For example, prior to integration with the other dis- 
ciplinary analyses, several modifications to CONSIZ were required 
so that accurate derivative information could be obtained in a rapid 
fashion. The lack of sensitivity information at design points other 
than the optimum is another drawback of the integrated approaches. 


Results and Discussion 

For each optimization method, solution of the single-stage-to- 
orbit problem is attempted from the two starting points listed in 
Table 1. The first starting point may be viewed as a good initial 
guess, whereas starting point 2 is quite poor. For each solution 
strategy, optimization was performed with the sequential quadratic 
programming algorithm NPSOL. 11 NPSOL uses a quasi-Newton 
method to approximate the Hessian of the Lagrangian. This algo- 
rithm is known to converge to local minimum for problems that are 
scaled properly and are twice continuously differentiable. 

Prior to optimization (with each solution strategy), the problem 
was scaled so that the design variables, constraints, and objective 
function were all of order one. However, in this investigation, no 
attempt was made to provide a twice continuously differentiable 
model. Sources that contribute to this nonsmoothness include the 
ascent pitch profile, which is modeled by discrete control points 
connected by linear segments, and the atmospheric properties (1976 
standard model). Hence, in this analysis, NPSOL is not expected to 
achieve a tightly converged, optimal solution, but rather to identify 
near-optimal solutions useful at the conceptual design level. All 
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computations were performed on a Silicon Graphics Challenge L 
machine outfitted with 100-MHz processors; CPU times quoted are 
based on use of this computational platform. 

Iterative-Loop Method 

The first approach toward an integrated multidisciplinary solu- 
tion strategy is depicted in Fig. 4. Here, an iterative loop is set up 
between the trajectory and weights and the sizing disciplines; val- 
ues of GLOW, S rcf , the base diameter, and the landed weight are 
used as loop convergence criteria. This type of strategy has been 
employed by other investigators 12 " 14 for solution of launch-vehicle 
design problems. Indeed, as an integrated design tool, this method 
should be viewed as the current way of doing business within the 
space transportation community. The popularity of the approach is 
largely a result of the tool-integration simplicity. Using the termi- 
nology of Ref. 15, this formulation may be referred to as “multi- 
disciplinary feasible” in that, for each set of design variables, the 
looped analyses return a design candidate that is consistent across 
disciplinary boundaries. 

Using this iterative-loop formulation, the single-stage-to-orbit 
launch-vehicle design problem was solved from each of the two 
starting points listed in Table 1. Each of these solutions is sum- 
marized in Table 2. Although the solution that began at point 2 
required 10 extra iterations, it is evident that both optimization pro- 
cesses converged to approximately the same point. At the solution, 
the vehicle liftoff thrust-to-weight ratio is the only variable on its 
bound (a liftoff T / W < 1 .2 is not allowed for abort reasons). 

The vehicle’s ascent profile is summarized in Figs. 5 and 6 for 
the solution that began at starting point 1. Flight begins with a 
400-ft vertical rise to clear the launch facility. This is followed by a 
maximum-pitch-rate segment in which the vehicle is trying to attain 
a maximum-lift orientation. The pitch rate is limited to 5 deg/s to re- 
flect control issues that are not modeled in this analysis. As shown in 
Fig. 6, during this segment of flight, the angle of attack ( a ) increases 
until it reaches the allowed maximum of 1 5 deg. Flying at this a , the 
normal force builds until a limit load is reached. The vehicle rides 
this normal-force boundary through peak dynamic-pressure, which 
for the optimum flight path is about 830 psf. Hence, the dynamic 
pressure limit of 1 000 psf is not active. During this phase of flight, at 
approximately 18 kft, the back-pressure losses are low enough that 
the dual-position nozzle is extended to gain propulsive efficiency. 
After the nozzle extension is completed, the vehicle acceleration ini- 
tially decreases as a result of flight through the transonic regime. At 



Fig, 4 Iterative-loop solution strategy. 


Table 2 Single-stage-to-orbit launch-vehicle design solution: 
summary for iterative-loop method 


From initial point 
1 2 


Looped variables 
GLOW, lb 

^rcf. ft 2 

Landed weight, lb 
Base diameter, ft 
Design variables 
Launch azimuth, deg 
Vehicle liftoff T/W 
% LH 2 , mode 1 
Mode 2 mixture ratio 
Nozzle area ratio, retracted 
Nozzle area ratio, extended 
Propulsion mode transition 
Mach number 

Nozzle transition Mach number 
Fz boundary duration, s 
Total trajectory time, s 
Performance 

Final objective (dry weight, lb) 
Number of iterations 
Number of analysis calls 
CPU time, h 


2. 1655e + 06 

2.1678e 4- 06 

4.0476e + 03 

4.0498e + 03 

2.2100e + 05 

2.2102e + 05 

27.916 

27.906 

42.458 

42.458 

1.2000 

1 .2000 

6.2587 

6.2190 

6.9986 

6.9985 

59.675 

59.687 

119.36 

119.30 

7.3137 

7.3608 

0.8904 

0.9011 

55.155 

60.248 

377.61 

377.62 

1 ,92676e + 05 

L92695e + 05 

65 

75 

10482 

13101 

13.3 

18.2 



Fig. 5 Optimal ascent altitude, velocity, and acceleration profiles ob- 
tained through application of the iterative-loop method: , altitude; 

, velocity; , acceleration. 


20 r 


15b 


10 - 

a, 

deg 

5 - 


Oh 


- 5 l 



Fig. 6 Optimal ascent pitch-angle, normal-force, and angle-of-attack 

profiles obtained through application of the iterative-loop method : , 

pitch angle; , F z \ , a. 


about 47 kft, as the dynamic pressure decreases, the vehicle comes 
off the normal-force boundary but continues to accelerate towards 
3 g. The vehicle reaches 3 g while in the dual-fuel propulsive mode 
and throttles down the engines to maintain this level of acceleration. 
In this analysis, transition of all seven engines from a dual-fuel to 
single-fuel mode is performed instantaneously. This results in the 
large decrease in acceleration shown in Fig. 5 at roughly 180 s and 
an altitude of about 164 kft. (A slightly better dry weight would 
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1 5 x 30-ft bay -i 
LH2tank^ \ 

A TnT^ 


LOX tank 






-181.3ft- 


Dry wt: 192 7 ktb 
Gross wt: 2. 16Mlb 


Fig. 7 Optimal single-stage-to-orbit vehicle description obtained 
through application of the iterative-loop method. 



a) Iterations 



Fig. 8 Objective function (dry weight) convergence history for the 
iterative-loop solution: a) first 10 iterations and b) final 65 iterations 
(note scale change). 

result if the engines were allowed to transition sequentially.) Oper- 
ating in a single-fuel mode (LH 2 ), the vehicle accelerates back to 
3 g and holds this acceleration until reaching orbit. 

Both solutions suggest that a vehicle dry weight of approximately 
192,700 lb is attainable for a 28-ft-diam configuration (181.3 ft in 
length) with a gross liftoff weight of about 2.165 x 10 6 lb and a 
wing area of 4050 ft 2 . Such a configuration is sketched in Fig. 7. 

Figures 8 and 9 show the objective function and nonlinear con- 
straint convergence histories beginning from the more demanding 
starting point. Note 1) the dry-weight scale change in going from 
Fig. 8a to Fig. 8b and 2) the semilog format of Fig. 9, which de- 
picts the norm of the constraint violations. From these figures, it 
is clear that after approximately 10 iterations the objective func- 
tion achieves a near-optimum value, whereas the constraints are not 
satisfied to the prescribed level (indicated by the nonlinear feasi- 
bility tolerance). These optimization characteristics are a result of 
the rather flat design space in which the most difficult task is to 
find a consistent vehicle model that satisfies all of the constraints. 
Lack of model smoothness is also evident in Fig. 9 and may have 
hampered the convergence rate. Note that at about iteration 60, after 
the optimizer has formed a good model of the design space, a large 



Nonlinear 

feasibility 

tolerance 


Fig. 9 Nonlinear-constraint convergence history for the iterative-loop 
solution. 


c 


OPTIMIZER 

Minimize J = dry weight 

x = (see Table 1) 

Subject to inflight and terminal constraints 


T rajectory 


Propulsion 


T 

Inflight & terminal 
constraints 






GLOW c Dry 

Sref c we, 9 ht 

Landed wt c 
Base diameter c 


Compatibility constraints: 

C(1) = 0 = GLOW c -GLOW 

C(2) = 0 = Sref c - Sref 

C(3) = 0 - Landed wt c - landed wt 

C(4) = 0 = Base diameter c - base diameter 


Fig. 10 Strategy for sequential compatibility-constraint solution. 


relative increase in dry weight occurs on achieving constraint fea- 
sibility. This feasibility is maintained as the dry weight is brought 
back down to the final value. 

In terms of the optimization strategy’s overall performance, one 
may note that a large number of analysis calls was required to achieve 
each of the solutions listed in Table 2. An analysis call is defined 
as one pass through the analysis loop, whereas an iteration includes 
all the calculations required to take a design step (loop convergence 
at the current design iterate, each point along the line-search, and 
all finite-differenced points). Given the fact that this problem was 
posed with 36 design variables and took approximately 70 iterations 
to reach the solution, one would expect slightly more than 2590 anal- 
ysis calls (for 70 objective functions and forward-differenced gradi- 
ents). However, these solutions required 4-5 times as many analysis 
calls. The source of this increased computational requirement is that 
while this strategy may provide the most straightforward approach 
to integrating an existing set of tools, multidisciplinary feasibility 
must be established at every point where either the objective func- 
tion or constraints are evaluated. In fact, this level of compatibility 
across disciplinary boundaries is only needed at the solution. In ad- 
dition to the added computational expense, this approach can lead to 
numerical difficulties when finite-differencing if the analysis loop 
is not converged tightly. In this investigation, a loop tolerance (E 
of Fig. 4) on the order of 1.0 x 10 12 was required to produce ac- 
curate derivatives. Without such a small tolerance, finite-difference 
errors led to spurious results. These computational difficulties are 
overcome in the next solution strategy. 
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Table 3 Single-stage-to-orbit launch-vehicle design solution: 
summary for sequential compatibility-constraint method 


From initial point 

i 2 

Design variables 


GLOW, lb 

2. 1675e + 06 

2. 1668e + 06 

■Vf. ft 2 

4.0500e + 03 

4.0497e + 03 

Landed weight, lb 

2.21 13e -P 05 

2.21 12e + 05 

Base diameter, ft 

27.928 

27.922 

Launch azimuth, deg 

42.328 

42.328 

Vehicle liftoff Tj W 

1 .2000 

1.2000 

% LH 2 , mode 1 

6.2614 

6.2596 

Mode 2 mixture ratio 

6.9893 

6.9996 

Nozzle area ratio, retracted 

59.541 

59.801 

Nozzle area ratio, extended 

1 19.08 

119.60 

Propulsion mode transition 



Mach number 

7.3062 

7.3117 

Nozzle transition Mach number 

0.8980 

0.9090 

Fz boundary duration, s 

63 400 

56.819 

Total trajectory time, s 

377.80 

377.60 

Performance 



Final objective (dry weight, lb) 

1 .92788e 4- 05 

1 ,92785e 4- 05 

Number of iterations 

70 

77 

Number of analysis calls 

3162 

3431 

CPU time, h 

4.0 

4.6 


Sequential Compatibility Constraint Method 

In the second approach, the iterative loop of Fig. 4 is replaced by 
the use of auxiliary variables and compatibility constraints. 16 ” 19 As 
shown in Fig. 1 0, an auxiliary variable and a compatibility constraint 
are added to the optimization-problem statement for each variable 
that is required as input to one discipline but is computed by another 
discipline later in the analysis sequence. Hence, 5 refl GLOW, the 
base diameter, and the landed weight are added as design variables 
in the Fig. 10 optimizer. In this manner, the iterative loop of Fig. 4 
is removed, and configuration control becomes an additional task 
of the optimizer. By satisfying these four compatibility constraints, 
a consistent vehicle model is guaranteed. However, as opposed to 
the approach of Fig. 4, compatibility is required at the solution only. 
This type of approach may be referred to as “simultaneous analysis 
and design,” since both a consistent and an optimum set of design 
variables are converged upon simultaneously. 20 

Using this formulation, solutions of the single-stage-to-orbit 
launch-vehicle design problem were obtained from each of the two 
starting points listed in Table 1. These solutions are summarized 
in Table 3. Once again, although the solution that began at starting 
point 2 required more iterations, both optimization paths converged 
to approximately the same point. Note that the solutions obtained 
through use of the compatibility-constraint approach are approx- 
imately 100 lb heavier than those found with the iterative-loop 
approach. This discrepancy results because the use of larger compat- 
ibility tolerances (1.0 x 10“ 5 ) was possible with this method. When 
these tolerances were reduced to 1.0 x 10~ 12 (the tolerance level 
that was required for the iterative-loop solutions), the compatibility- 
constraint solutions were within 20 lb of the looped solutions. Recall 
that in the iterative-loop approach, a compatibility tolerance on the 
order of 1 0“ 5 was not possible because of the loss of finite-difference 
accuracy. Hence, a major advantage of the compatibility-constraint 
approach is that accurate derivatives may be obtained without re- 
quiring strict convergence of the disciplinary models. This flexibil- 
ity may be used to realize reduced computational requirements in 
preliminary design analyses where coarser approximations are ade- 
quate, while maintaining the capability for more accurate solutions 
as the design progresses. Another computational advantage results 
from the fact that a consistent vehicle model is only required at the 
solution. Therefore, disciplinary-model convergence is not required 
at each design iterate, each point along the line search, and each 
finite-difference point. 

These computational advantages are clear in the number of 
analysis calls (and corresponding CPU time) required to obtain 
each solution. Compared to the iterative-loop method, the se- 
quential compatibility-constraint approach requires only 25% of 


the computational time to reach approximately the same solution 
(within 0.05% in dry weight). Recall that in this approach we have 
40 design variables (the original 36, plus the four auxiliary variables 
used to establish compatibility). Hence, for 70 iterations, one would 
expect slightly more than 2870 analysis calls. This estimate matches 
well with the actual number of analysis calls listed in Table 3. The 
only disadvantage to this approach relative to the looped method 
would be in a situation where solution convergence is not achieved. 
In that case, one could end up with no valid design information, since 
model convergence is not guaranteed across disciplinary boundaries. 
Note that this situation, which could arise in a poorly scaled prob- 
lem, may be mitigated by adjusting the scaling strategy 21 or through 
application of a feasible-point optimization method. 

The ascent-trajectory time histories for this case were nearly iden- 
tical to those obtained with the previous solution strategy (Figs. 5 
and 6). Convergence histories are shown for the objective function 
and norm of the nonlinear constraints beginning from starting point 
2 in Figs. 11 and 12. Note the dry-weight scale change in going 
from Fig. 1 la to Fig. 1 lb, and the semilog format of Fig. 12. From 
these figures, it is clear that after approximately five iterations the 
objective function achieves a near-optimum value, and after ap- 
proximately 20 iterations the compatibility constraints are satisfied. 
However, as in the looped approach, it is the convergence of the 
highly nonlinear trajectory constraints that is causing computational 
difficulty. Once again, the lack of model smoothness is evident in 
Fig. 12 and may have hampered the convergence rate. Comparing 
Figs. 1 1 and 12 with Figs. 8 and 9, it is clear that the use of compat- 
ibility constraints provides another advantage — a smoother design 
space. That is, by eliminating the loop convergence requirement, a 
smoother analysis results. This same effect was noted for conver- 
gence from starting point 1 . 

Note that the approach of Fig. 10 still requires a sequential so- 
lution. That is, the propulsion discipline is analyzed first, followed 
by the trajectory discipline, and then finally the weight and siz- 
ing calculations, A sequential solution is required because only the 
feedback portions of the iterative loop were broken; feedforward 




b) 

Iterations 

Fig. 11 Objective-function (dry weight) convergence history for the 
sequential compatibility-constraint solution: a) first 10 iterations and 
b) final 67 iterations (note scale change). 
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Fig. 12 Nonlinear-constraint convergence history for the sequential 
compatibility-constraint solution. 

from one discipline to the next is still performed (e.g., the mass 
ratio is computed by the trajectory discipline and used as an in- 
put for the weight and sizing analysis). This differs from the ap- 
proach suggested in Ref. 22, where removal of both the feedforward 
and feedback communication requirements yields a parallel analysis 
strategy. A sequential analysis strategy was selected for solution of 
this problem because it required the addition of only tour auxiliary 
variables and compatibility constraints. In comparison, a parallel 
strategy would have required an additional 21 variables and con- 
straints, nearly doubling the size of the original problem. Although 
increasing the problem size does not always result in an increased 
solution time (the optimization algorithm, the order and sparseness 
of the design space, and the implementation strategy are all major 
factors), 23,24 for a dense, sequential quadratic programming algo- 
rithm (such as NPSOL) the solution time is roughly proportional to 
the problem size. 

Summary 

The present investigation focused on the development of a rapid 
multidisciplinary analysis and optimization capability for launch- 
vehicle design. In particular, two multidisciplinary optimization 
strategies were implemented and evaluated for solution of a singie- 
stage-to-orbit launch-vehicle design problem. Weights and sizing, 
propulsion, and trajectory issues were directly addressed, and the 
need to maintain a consistent vehicle model across the disciplines 
was discussed. The approaches were referred to as the iterative-loop 
method and the sequential compatibility-constraint approach. 

Both solution strategies require some disciplinary tool integra- 
tion and were shown to reach approximately the same solution from 
two different starting points. Each of these solutions suggests that a 
dual-fuel, single-stage-to-orbit vehicle with a dry weight of approx- 
imately 1.927 x 10 5 lb, gross liftoff weight of 2. 165 x 10 6 lb, and 
length of 181 ft is attainable. 

Currently, within the launch- vehicle design community, methods 
of the general form of the iterative-loop approach are predominantly 
used. In the present investigation, the manner in which disciplinary 
coupling is treated was shown to have a direct effect on optimiza- 
tion convergence and the computational effort required to reach 
the solution. This analysis has shown that use of the sequential 
compatibility-constraint approach has several advantages relative to 
the iterative -loop approach. These advantages include: 1 ) being 3-4 
times more computationally efficient (for this problem), 2) provid- 
ing greater flexibility in the way in which consistency is maintained 
across disciplinary boundaries, and 3) a smoother design space. 
The only disadvantage of the compatibility-constraint approach is 
in situations when the optimizer terminates without reaching the 
solution on account of poor scaling or model nonsmoothness. Be- 
cause multidisciplinary feasibility is only guaranteed at a solution 
in this approach, the design information could be invalid. How- 
ever, the likelihood of optimization termination without constraint 


satisfaction could be lessened with the use of a more robust scaling 
strategy or a feasible-point optimization algorithm. 
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